ON UPPER LIMITS FOR GRAVITATIONAL RADIATION 



P.Astone^ and G. Pizzella^ 
^ INFN, Sezione di Roma 
^ University of Rome Tor Vergata and INFN, Laboratori Nazionali di Frascati 
. P.O. Box 13, 1-00044 Frascati, Italy 

O 

o 

(N 

a 

(N 

(N 
> 

2 ■ Abstract 

. A procedure with a Bayesan approach for calculating upper limits to gravitational 

I wave bursts from coincidence experiments with multiple detectors is described. 
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1 Introduction 

After the initial experiments with room temperature resonant detectors, the new 
generation of cryogenic gravitational wave (GW) antennas entered long term data taking 
operation in 1990 (EXPLORER 0), in 1991 (ALLEGRO 0), in 1993 (NIOBE 0), in 
1994 (NAUTILUS i) and in 1997 (AURIGA §). 

Searches for coincident events between detectors have been performed. Between 
EXPLORER and NAUTILUS and between EXPLORER and NIOBE in the years 1995 
and 1996 §. Between ALLEGRO and EXPLORER with data recorded in 1991 [0]. In 
both cases no significative coincidence excesses were found and an upper limit to GW 
bursts was calculated 0. 

However, the upper limit determination has been done under the hidden hypothesis 
that the signal-to-noise ratio (SNR) is very large. According to theoretical estimations the 
signals expected from cosmic GW sources are extremely feeble, so small that extremely 
sensitive detectors are needed. In fact, according to present knowledge, the detectors 
available today have not yet reached the sensitivity to detect even a few events per year. 

Thus it is important to study the problem of the upper limit determination in the 
cases the SNRs of the observed events are not large. In order to do this we have to discuss 
our definition of event. 

The raw data from a resonant GW detector are filtered with a filter matched to 
short bursts We describe now in more detail the procedure used for the GW detectors 
of the Rome group, EXPLORER and NAUTILUS. 

After the filtering of the raw-data, events are extracted as follows. Be x{t) the fil- 
tered output of the electromechanical transducer which converts the mechanical vibrations 
of the bar in electrical signals. This quantity is normalized, using the detector calibration, 
such that its square gives the energy innovation Ej of the oscillation for each sample, 
expressed in kelvin units. In absence of signals, for well behaved noise due only to the 
thermal motion of the bar and to the electronic noise of the amplifier, the distribution of 
x{t) is normal with zero mean. The variance (average value of the square of x{t)) is called 
effective temperature and is indicated with T^ff. The distribution of x{t) is 



For extracting events (in absence of signals the events are just due to noise) we set a 
threshold in terms of a critical ratio defined by 
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The threshold is set at a value CR such to obtain, in presence of thermal and 
electronic noise alone, a number of events which can be easily exchanged among the other 
groups who participate to the data exchange. For about one hundred events per day the 
threshold corresponds to an energy Et = 19.5 T^jf. 

We calculate now the theoretical probability to detect a signal with a given SNR, 
in presence of a well behaved Gaussian noise. We put y = {s + xY where s = V SNR is 
the signal we look for and x is the gaussian noise. We obtain easily [§] 

nOO -| SNR 

probability's NR) = / -^=6'^^^^ cosh{Jy ■ SNR)dy (4) 

J SNRt V27r2/ 

We put SNRt = 19.5 for the present EXPLORER and NAUTILUS detectors. 

2 Upper limit determination 

We consider M detectors and search for M-fold coincidences over a total period of 
time tm during which all detectors are in operation. Be n the average number of accidental 
coincidences (due to chance) and Uc the number of coincidences which are found within 
a given time window. 

For events which have a Poissonian distribution in time the expected average number 
of M-fold accidental coincidences is given by 



1,M 

T, r A,f_1 

n 
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where is the event density of the k^^ detector. 

The accidental coincidence distribution can be estimated experimentally by proper 



shifting []IU| the event occurrence times of each detector. In the case of Poissonian dis- 
tribution the average number of the M-fold accidental coincidences coincides with that 
given by eq. ^. The comparison between ric and n allows to reach some conclusion about 
the detection of GW or to establish an upper limit to their existence. 

In paper [0 and in the previous paper |Tl| the upper limit has been estimated 



as follows. It has been found that, for various energy levels of the observed events, the 
number ric was smaller than or did not exceeded significantly n. Such numbers nc, one 
for each energy level, were used for calculating the upper limit. A Poissonian distribution 
of the number of the observed events was considered together with the hypothesis of 
an isotropic distribution in the sky of the GW sources. The value of h (adimensional 
perturbation of the metric tensor) was then derived from the energy levels, using the 
detector cross-section for gravitational waves. 

This procedure can be objected on two points: 



a)The most important point is that, as shown in |T^, for SNR small and up to values of a 
few dozens, the energy of an event is not the energy of the GW absorbed by the detector. 
This means that we cannot deduce the value of h directly from the energy levels of the 
observed events; 
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b)In addition, the efficiency of detection, again for SNR values up to one or two dozens, 
is rather smaller than unity, and this changes the upper limit, particularly at small SNR. 

We introduce a new procedure for estimating the upper limit, which circumvents 
the difficulties indicated in the above two points. 

The problem to determine the upper limit has been discussed in several papers. In 
particular in paper [|15|, as indicated by the PDG, [l^ and, more recently, in paper fl% . 



According to [0 the upper hmit can be calculated using the relative belief updating ratio 



R{nGw,n,,n) = e-"«-(l + ^)"^ (6) 

n 

referring to a given period of data taking. This function is proportional to the likelihood 
and it allows to infer the probability to have now signals for given priors (using the Bayes's 
theorem). It has already been used in High Energy Physics 0, [l^ . 
We calculate the upper limit by solving the equation 

R{nGw,nc,n) = 0.05 (7) 

We remark that 5% does not represent a probability but it is an useful way to put a limit 
independently on the priorsj^. 

Eq. 1^ has a very interesting solution. Putting Uc = we find uqw = 2.99, inde- 



pendent on the value of the background n. If we use the calculations of ref. |15| we find 
that, for ric = and n = 0, the upper limit is 3.09 (almost identical to the previous one) 
but it decreases for increasing n. The reason for this different behavior is due to the non- 



Bayesan character of the calculations made in |TB[, as we discuss in the following. 

Suppose we have ric = and n ^ 0. This certainly means that the number of acci- 
dentals, whose average value can be determined with any desired accuracy, has undergone 
a fluctuation. For larger n values, smaller is the (a priori) probability that such fluctua- 
tion occur. Thus one could reason that it is less likely that a number new be associated 
to a large value of n, since the observation gave Uc = 0. 



According to the Bayesan approach instead, as discussed in [0, one cannot ignore 
the fact that the observation ric = had already being made at the time the estimation of 
the upper limit is considered. The Bayesan approach requires that, given Uc = and n 7^ 0, 
one evaluate the chance that a number new of signals exist. This chance of a possible 
signal is referred to the observation already made and, rather obviously, it cannot depend 
on the previous fluctuation of the background, since the presence of a signal cannot be 
related to the background due to the detector. Mathematically, it is easy to demonstrate, 
using the results obtained in that due to the Poissonian character of the number of 
accidentals this relative chance (for nc = 0) is indeed independent on n. 

It can be seen, comparing the results of [|T5| with those of [0, that the Bayesan 
upper limits are for all values of tt-c and n (except nc = n = 0) , greater than those obtained 



To avoid confusion we shall continue to use the words upper limit, although it would be more 
appropriate to call it standard sensitivity bound |20| . 
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Figure 1: Number of GW signals expected for the sensitivity limit of 5% versus the 
number of coincidences equal to the average number of accidentals. 

with the non-Bayesan procedure. In our opinion the Bayesan approach has to be preferred, 
and so we do in this paper. 

If we have ric ^ then we apply eq. |^. It is interesting to show the result for the 
case Uf. = n ^ for the standard sensitivity bound of 5%. The result is given in fig.|I| We 
note that for ric = n and new << ^ eq. ^ can be approximated with 

now ~ (8) 

From the result shown in fig.|l] it appears evident that the lowest upper limit is 
obtained for tt-c ~ ti ~ 0. In order to obtain fi ~ one can raise the threshold used for 
determining the events. However in doing this one diminish the efficiency of detection, 
as shown in eq.^. Whether the procedure to raise the threshold is convenient or not, it 
depends on the numerical effects of the two competing operations. Certainly for large GW 
signals, when the detection efficiency is always unity, it is much better to have a threshold 
that gives n = 0. For smaller signals one has to consider specific cases. However it can 
be seen that in the most interesting cases it is better to raise the threshold until we get 
n ~ 0. This will be shown in the section where we reconsider the upper limit obtained 
with ALLEGRO and EXPLORER in 1991 §. 

In the estimation of the upper limit we consider the efficiency of detection, which 
we indicate with ek{SNR) where k refers to the k^^ detector. For EXPLORER and NAU- 
TILUS the theoretical efficiency is obtained from eq. ^. 

We must relate the h values of the GW to the energy E absorbed by the detectors. 
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We have to consider that the absorbed energy depends on the direction of the impinging 
GW and on its polarization. For taking care of the various polarization we use the average 
value dividing the cross section by a factor of two. We then have 



h = 1.13 10^'VS (9) 

with the energy E expressed in kelvin unit. This formula is valid only if the GW arrives 
perpendicularly to the detector axis {9 = 90°). For a given direction we calculate the 
absorbed energy using the sin{9)^ dependency. We also consider that for an isotropic 
distribution of sources the number of possible GW impinging directions is proportional 
to siniOy. 

The procedure for calculating the upper limit is accomplished thru the following 
points: 

a) consider various values of /i; 

b) assume an isotropic distribution of the GW sources; 

c) for each direction 9 and for each h calculate the absorbed energy E{9) by means of eq. 
^ and the sin^{9) dependency; 

d) for each detector calculate the SNR for the adsorbed energy by taking into consideration 
the noise T^jf^k- 

SNR,{9) = ^, k = l,..,M (10) 

e) using the individual efficiencies ek{SNRk{9)) consider the total efficiency et{9) = 
U]:''e,{SNR,{9)); 

f) integrate et{9) over 9 with the weight sin'^{9), because of the assumed isotropic distri- 
bution of the sources; 

g) from eq.p, given ric and n, we obtain new- We then divide new by the result of point 
f) and obtain for each value of h the upper limit during the measuring time tm- 

We remark that in this case we have not used the energy of the observed events, as 
done instead previously |[Tl|, 0]. 

The total efficiency is calculated with the following eq. |Tl]. 

^^^^(^^ ^ j}Ul:''ek{SNR,{9))stn\9)d9 ^^^^ 

4 

For more clarity we show in table |l] some of the steps needed for our calculation, 
using two parallel detectors and Uc = 0. We use the efficiency given by eq. |[ valid for a 
well behaved nois^^ 

3 Ricalculation of the upper limit with the data of ALLEGRO and 
EXPLORER in 1991 

In a previous paper ||^ the upper limit for GW bursts was calculated, using the data 
recorded by ALLEGRO and EXPLORER in 1991. We wish now to recalculate the upper 
limit according the considerations discussed in this paper. 
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Table 1: Procedure for calculating the upper limit with two detectors. We assume that 
one detector has noise Tejj = 1 mK, the other one has noise Te/j = 2 mK. For each 
value of h we give: maximum energy adsorbed by the detector (for sin'^{6) = 1), SNR 
and efficiency of detection for each detector, total weighted efficiency (having considered 
an isotropic distribution of the GW sources. Due to the angular weighting etotai < ^a^b)- 
The upper limit is given by 

^total 



h 


Eabs 


Detector A 




Detector B 






upper limit 




[mK] 


SNRa 




SNRb 




^tot 
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31 


31 


0.88 
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0.32 


0.12 


16 
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70 


70 
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35 


0.93 


0.56 


4.3 


4 


126 


126 
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63 


1 


0.76 


3.6 


10 


783 


783 


1 


392 


1 


0.95 


3.1 



In 1991 the EXPLORER data filtering was done differently from that described in 
this Introduction. For both ALLEGRO and EXPLORER the output of the electrome- 
chanical transducer was sent to lock-ins referred to the frequencies of the resonant modes. 
Then the outputs of the lock-ins (in phase and in quadrature) were filtered searching for 
delta-like signals and combined for obtaining the energy innovation, which we still indi- 
cate with Ef. In this case the probability to have an event (above threshold SNRt) due 



Here y = t^t^, Iq is the modified Bessel function of order zero, and the noise temperature 



to a signal with given SNR is obtained (see ref. P, [T^) with the following equation: 

poo 

probability{SNR) = / g-(5ivi?+j;) j^^2^y ■ SNR)dy (12) 

J SNRt 

iL 

•ff 

Tf,ff is the average value of the energy innovation Ef. 

We recall that in a time period of 123 days 70 coincidences were found with a 
background of 59.3. For extracting the events the ALLEGRO threshold was SNRt = 11.5 
with a noise temperature T^ff ~ 8 mK. For EXPLORER the threshold was SNRt = 10 
also with Teff ~ 8 mK. Applying eq.^ we find an upper limit of now = 37 over the 123 
days. 

According to the previous considerations we can raise the event threshold, say for 
EXPLORER, in order to reduce the number of accidentals. For instance, for a threshold 
SNRt = 24 we get ric = I and n = 0.74, obtaining, from eq. |^, the value uqw = 4.8. 

Thus the procedure for calculating the upper limit with the Bayesan approach when 
we have data at various thresholds, including cases with Uc and n different from zero, is 
the following. 

Start with ric and n for various thresholds and use eq.^ for obtaining new at each 
threshold. Calculate the upper limit for various values of h as shown in the previous 
section. For each h take as upper limit the smallest value among those obtained by varying 

The real data often show a non gaussian behaviour. In this case the efficiency differs from the theo- 
retical one given by eq.^, but one can easily make use of the efficiency experimentally measured. 
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Figure 2: The asterisks indicate the upper hmit calculated in |^. The other line indicates 
the upper limit evaluated with the Bayesan approach. 

the threshold. Clearly at large h values, when we get Uc = 0, the upper limit is, for the 
entire period of time, now = 2.99. 

The result is shown in fig.^] together with that obtained previously in . It turns 
out that the two upper limits are similar. 

The reason for this is due to the fact that in applying the previous algorithm we 
started from an energy level higher than the largest energy of the detected (accidental) 
coincidences, thus obtaining, at this level [ric = fi = 0) an upper limit of 3.09 very close 
to the value 2.99 obtained with the Bayesan approach. The similarity of the results at 
lower h values is accidental. In the previous algorithm the increase at lower h is due only 
to the increase of the number n of accidentals. In the present algorithm the increase is 
due to the smaller efficiency of detection and to the increase in uqw which roughly goes 
with (eq-|)- 

In spite of the similar numerical results, we believe that the procedure proposed here 
which does not extract the value of h from the energy levels of the accidental coincidences 
and it uses the Bayesan approach is methodologically more correct. 

4 Discussion 

The best upper limit which can be obtained with an array of M identical parallel 
detectors in M^' coincidence cannot go below the value 2.99, because this is the upper 
limit WTf when one finds zero coincidences independently on the background. 



The basic advantage in using many detectors comes from the fact that with many 
detectors it is easier to obtain ~ 0, and thus (in absence of GW) = 0. Because of 
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the Poisson distributions, the average number of accidental coincidences for M detectors 
in a time window ±w is given by eq.^ On the time scale of 1 second (w=l s) it turns out 
that Uk << 1. By increasing the number of detectors one obtains smaller values of n, thus 
approaching the requirement to have ric = and then the lowest possible upper limit. 

This is certainly true at large h values, where the detection efficiency for all detectors 
is unity. The result, as shown in fig.^], is a plateau. Instead it might be convenient at low 
h values to use the two most sensitive detectors, in order to have the largest possible 
efficiency of detection. The overall upper limit is then obtained by taking the smallest 
ones among the values of the various upper limit determinations. 

The above procedure can be easily adjusted to the more general case of any distri- 
bution of the GW sources, and of non-parallel detectors. 
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